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Four-parameter Beta Compound Binomial 
Abstract 

This paper presents a detailed derivation of method of moments estimate for the four- 
parameter beta compound binomial strong tnie score model. A procedure is presented to 
deal with the case in which the usual method of moments estimates do not exist or result 
in invalid parameter estimates. The results presented regarding the method of moments 
estimates are used to derive formulas for computing classification consistency indices imder 
the four-parameter beta compound binomial model. 

Acknowledgement. The author thanks Robert L. Breraian for carefully reading an earlier 
version of this paper and offering helpful suggestions and comments. 
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Four- parameter Beta Compomid Binomial 
The first part of this paper discusses estimation of the four-parameter beta compound 
binomial model by the method of moments. Most of the material in the first part of 
this paper is a restatement of material presented in Lord (1964, 1965) with some added 
detaik The second part of this paper uses the results presented in the first part to obtain 
formulas for computing the classification consistency indexes described by Hanson and 
Brennan (1990). 

The purpose of this paper is to provide a detailed description of the procedures used in 
computer programs written by the author which compute estimates for the four-parameter 
beta compound binomial model and indexes of classification consistency based on the four- 
parameter beta compound binomial model. 

Estimation of the Four-parameter Beta Compound Binomial Model 

It is assumed that the test score to be modeled is the sum of K dichotomously scored 
test items (tliis score is referred to as the raw or observed test score). The probabiHty 
that the raw score random variable X in the population of interest equals i {i = 0, . . . , A'), 
under the four-parameter beta compound binomial model is given by 



(1) 



where r is the proportion correct true score. (For simphcity of notation the dependence 
of the marginal test score distribution {Pr(X - i)] on the parameters k, a, ^, /, and u 
is not denoted.) The true score distribution [^(r j a,^J,u)] is assumed to belong to the 
four-parameter beta family of distributions. The four-parameter beta distribution is a 
generalization of the usual beta distribution that in addition to the two shape parameters 
(a > 0 and l3 > 0) has parameters for the lower {/) and upper (h) Hmits of the distribution 
(0 < / < tx < 1). The four-parameter beta density function [defined on the interval (/,«)] 
is 

g{r I cijij.u) 



(t/-/)«+^-iS(o,^) ' (2) 
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Eour-pbrametcr Beta Compound Binomial 
where B(a, 0) is the beta function, which is related to the gamma function \V{x)] by 

The conditional error distribution [Pr(A" = i | T,k)] is assumed to be Lord's (1965) 
two term approximation to the compound binomial distribution. The probability density 
function of the two-term approximation to the compound binomial distribution is (Lord, 
1965, Equation 5) 

VTiX = i\r,k)=pii\K,T) 

-Ml - r)[pii. i Jv - 2,r) - 2p{i - 1 | A' - 2,r) ~2\h'~ 2,r)] , (3) 

where p{i | n, r) = Pr(Z = i | n, r ) where the distribution of Z is binomial with parameters 
n (the number of binomial events) and r (the binomial probability). 

The two-term approximation to the compound binomial distribution given by Equa- 
tion 3 involves the parameter k {note that this is distinct from upper case A" which is 
used in this paper tc donate the muiiber of items on the test), in addition to the binoniial 
parameter. When k = 0 the conditional distribution of X given r is binomial. Lord (1965) 
gives a method of estimating k by setting the theoretical value of the average error variance 
(or reliability) under the two term approximation of the compound binomial distribution 
(which is a function of k) equal to an estimate of the average error variance (or reliabil- 
ity). The number-correct true score \'ariance, assinning the condititwial error distribution 
is Lord's two-term approximation to the compound binomial distribution, is (Lord, 19G5. 
Equation 39) 

crl-{K-2k)nAK-f'r) . 
A'(A'-1) + 2A: 

where fii and are the raw sc<)re mean and variance. Subtracting Efiuation 4 from aj 
gives the average error variance (denoted cr^) under the two term approximation of the 
c^nnpoimd binomiid distribution 

, aj [K{K -l) + 2k- K-'] -h (A - 2A')//,( A - ft,) 

A'(A'-1) + 2A- • 



Four-parameter Beta Compound Binomial 

Solving Equation 5 for k gives 

2[M.(A:-^,)-(<,|-.r?)] • ^' 

Given an estimate of the awrage error variance, <tJ, and estimates of fix and <tJ in Equation 
6, an estimate of k can be calculated. The value of Jb as a function of the reliability (/?), raw 
score mean, and raw score variance is given by substituting iT^(l — p) for tr^ in Equation 
6. 

The conditional error distribution given by Equation 3 may not be a probability 
distribution since for some values of r and i it is possible that Pt{X = t | r, jt) < 0. Lord 
(1965) states that for usual values of A' and k negative probabilities are typically negligible 
for .01 < r < .99 so that for practical purposes it is appropriate to treat Pr(A' = i ) r, k) 
given in Equation 3 as a probability distribution. 

After a value of k has been determined there are two steps involved in the estimation of 
the observed score distribution under the four-parameter beta compound binomial model 
using the method of moments. First, the parameters of the four-parameter beta true score 
distribution are estimated, and second these parameter estimates are used to produce the 
estimated observed score distribution. 

Estimation of IVue Score Distribution 

Under the assumption that the conditional error distribution is given by the two term 
api)roximation to the compound binomial distribution. Lord {19G5) produces a formula 
which gives the non-central moments (moments about zero) of the proportion correct true 
score distribution in terms of the factorial moments of the observed score distribution. The 
first moment of the proportion correct true score distribution can be written as 

1 

fhT = ^P'm., (7) 

where /i,^ is the i-th central uionient of the proportion correct true score distribution and 
/'li)x factorial moment of the niunher correct observed score tiistribution (the 




Four-paraxneter BvtR Coinpoimd Bmomiid 

i-ih factorial moment of a random variable X is defined as EjA'C-Y -f 1) • • - (-^ + ^ ~ 1)])- 
The i-th (i > 2) non-central moment of the proportion corrrel t' ue score distribution 
is (Lord, 1965, Equation 37) 

A'(A' - 1) + ibiW ' ^ ' 

where for integers j and /, j^"^ = i(i - 1) • • • 0' - ' + !)• Given estimates of the first four 
factorial moments of the observed score distribution {these can be obtained from the central 
or non-central observed score moments, see Kendall & Stuart, 1977, page 66), estimates of 
the mean and second through fourth non-central moments of the proportion correct true 
score distribution can be obtained from Equations 7 and 8. 

The estimates of the first four true score moments are used to produce method of 
moments estimates of the parameters of the four-parameter beta true score distribution. 
The four parameters to be estimated are the two shape parameters (a and /?), the lower 
limit of the distribution (/) and the upper limit of the distribution (u). If the i>n>i>ortion 
correct true score distribution is a member of the foiur-parameter beta family then its 
mean, variance, skewness and kurtosis are given by (Johnson L Kotz, 1970, pages 40-44) 



^2r = (U - /) : 



(9) 



■y - -2(3 a) V^^-^^+^ 

where is the 2-th central moment (moment about the mean) of the proportion correct 
true score distribution. The central moments of Equation 9 can be written in tenns of tli<' 
non-central moments of Equation 8 (see Kendall k Stuart, 1977, page 58). 
Let 



(6 -f 3-71-274) 

9 



r= ^ urn 
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Four-parajnetcr Beta Compound Binominl 
Substituting the values for 73 and 74 given in Equation 9 into Equation 10 and simplifying 
gives r = + Substituting r for a + in the expression for 7^ given in Equation 9 and 
solving for q0 gives 

6r^ri-l) 

{r + 2)(r + 3)74-3(r-6)(r-|-l)- ^^^^ 
Solving for a using the expressions a + /? = r and Equation 11 gives the two solutions 

2V V (»' + 2)(r + 3)7,-3{r-6)(r + l) y * ^^2) 

If the st>lutions given in Equation 12 are real then one of these solutions will be the value 
of a and the other solution will be the value of /?. If the skewness (73 ) is positive then the 
larger solution will be the value of 0, otherwise the larger solution wiU be the value of a. 

Substituting estimated moments for population moments in Equation 12 gives method 
of moments estimates of the parameters « and 0. Solving the expressions for and /i2r 
in Equation 9 for / and u gives 

' - /^Ir - ■ ■7==y 

yap 

Substituting estimates of thr and the method of moments estimates of a and 0 from 
Equation 12 in Equatif«i 13 gives mrthod of moments estimates of / and u. 

The mrthod of moments estimates of the four parameters do not exist when the 
solutions for o and 0 given in Equation 12 art> not real. Even if solutions exi«t the estimates 
of some parameters may not be valid (i.e., « > 1, / < 0, a < 0, 0 < 0). When the method 
of moments solution using the first four nK,mrnts does not exist or onv or more of the 
estimated parameters k invnlirl it is .suggested an estimation procedure l)e usetl in which 
the first tliree moments nre fit to determin<. three of the four paramett^s {a, 0 and u) and 
the remaining ]>aramet«T (/) is chosen such that tli<. kurtosis of the fitted distribution is as 
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Four-paranieter Beta Compoxmd Binomial 

near as possible to kurtosis as calculated directly from the observed score moments (using 
Equation 8). To implement this procedure, method of moments estimates of a, and u 
(given a specijSed ^'alue of /) based on the first three moments of the proportion correct 
true score distribution are needed. 

To derive method of moments estimates of a, /9 and u (given a specified value of /) 
the following definitions are made 

6= t = 1,2,3. (14) 

Using Equation 14^ ^3 can be written as a function of ^1 and 

The first two non-central moments of the proportion correct true score distribution can be 
written as 

(16) 

/4, = /2 + 2/(«-/)6 + («-0'6 6. 

Solving Equation 16 for 6 ^"^^ ^2 gives 

(/,^ -/)(«_/) • 

Substituting the expressions for and given in Equation 17 into the expR-ysion for ^3 
given in Equation 15 and simplifying yields 

The third non-central monieut of the proportion correct tnie score <ljslribntion is 

= + 3f'{u /)6 -f 3l{u - If 66 -f {'^ - /)' ^iUa . (10) 

6 



Four-parameter Bctn Compound BinoiniaJ 

Substituting the expressions for ^ and 6 from Equation 17 into Equation 19 and simpli- 
fying produces 

ri, = + 3/2(;.;, - /) + 3/[4, -F- 2l{ti\^ - /)] + [^i^ - /2 _ 2/(,.; , - /)] in - /)6 . (20) 

The only expression involving u in Equation 20 is (u - Substituting the right hand 
side of Equation 18 for (u - /)5 in Equation 20, solving for u and simplifying gives 

Substituting estimated moments in Equation 21 in place of the population moments gives 
a method of moments estimate of u (/ is assumed to be specifii-d). 

Let r* be a transformation of the proportion correct true score given by 



• III 



T' = 



(22) 



Since r has a four- parameter beta distribution, r* has a two-parameter beta distribution 
with parameter a and /I The parameters a and ^ are given in terms of the non-central 
moments of tlie distribution of r* as 

/'2r--(M..)' 

a = (i-M.0(/4.>-/4.O ^23) 

Given vcihies of u and /, the non-central moments of the distribution of r* can be obtained 
from the non-central moments of the distribution of r. Substituting estimated moments 
of the distribtition of r' m Equation 23 in place of the population moments gives method 
of moments estimates of o and 

Given a vahie of /, Equations 21 and 23 can be used to obtain method of moments 
estimates of w. n. and If the method of moments estimates of /, a and ^ given by 
Equations 12 and 13 <lo not exist or are not %-jdi(i, then a solution is selected such that 
the first tiiree nionients are fit (which determines u, an<i li) and / is chosen such that 
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Four-parameter Beta Compound Binomial 

the fitted true score kurtosis is as close as possible to the true score kurtosis as calculated 
directly from the observed score moments (using Equation 8). 

The expressions for a and /? given in E^u^cion 23 can be written in terms of / and the 
first three non-central moments of the proportion correct true score distribution (/fj^, /ijr 
and These values of q and ;9 can be used in the expression for the kurtosis given in 
Equation 9 to compute the fitted true score kurtosis as a function of /. The function of / 
given by the squared difference of fitted true score kurtosis (a function of /) and the true 
score kurtosis calculated directly fr-.m the observed score moments (tising Equation 8) is 
to be mininuzed. The minimization is constrained to be over those values of / for which 
valid method of moments estimates of «, a and 0 using Equations 21 and 23 exist. 

A method of finding the value of / that solves this constrained optimization problem is 
to first compute the two solutions in which / and u are on the boundary of the parameter 
space (the solution for which / = 0 and the solution for which u ~ I). Either or both of 
these solutions may not be a valid solution. The solution for which / = 0 can be computed 
using Equations 21 and 23 . The solution for u = 1 can be computed by solving Equation 
21 for / giving 

Substituting u = 1 in Equation 24 gives a method of moments estimate of /. Tliis estimate 
of \ along with m = 1 are used in Equation 23 to give estimates of a and /3. For each 
of these two solutions (the solution for which I = 0 and the solution for which w = 1) 
the fitted kurtosis is calculated (assuming the solutions arc valid). The solution with the 
smallest squared difference in the fitted kurtosis and the kurtosis calculated using Etjuation 
8 (this squared difference will be referred to as the squared kurtosis difference) is used as 
the initial solution. A grid search (Thisted, 1988, page 200) is then used for values of / > 0 
to search for a solution with a smaller squared kurtosis difference than the initial solution. 
For almost all the situations that have been encountered in practice rilher the sohition 
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Four-parameter Beta Compound Binomial 

for which / = 0 or the solution for which « = 1 produces the smallest squared kurtosis 
difference. There are examples, though, in which a solution with / > 0 and u < I can 
give a smaller squared kurtosis difference than either the solution for which / = 0 or the 
solution for wliich u = 1. 

Estimation of Observed Score Distribution 

Tliis section discusses estimation of the observed score distribution under the four- 
parameter beta compound binomial model assuming a value of k has been determined and 
estimates of a, /, and u have been calculated. The derivations presented in this section 
closely follow those presented in Lord (1964) with some added details. These derivations 
and the resulting formula for computing the oteerved score probabilities are not presented 
in Lord (1965). 

Casel: ^ 0. The case is first considered in which the conditional error distribution 
[Pr(X ^ i I r,k)] in Equation 3 is binomial (A: = 0). In this case the observed score 
distribution is 



Pr(.V = i) = J' (f )..(! - .) 



K-i 



dr 



(25) 



Let T* be defined as in Equation 22 so that r = (« - /)r* + / and dr = {u~ l)dr*, MaJcing 
the change of variable from r to r* in Equation 25 gives 



Pr( 



{u-l)dT* 



(26) 



Simplifying tlic expression of the true score density times u - / in Equation 26 and substi- 
tuting (1 - u) + (« - /) (1 _ r*) for l-{u-l)r' - I gives 

•(r-)"-^(l-r*)^-^/r\ (27) 
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Four-parameter Beta Compound Binomial 
Using the binomial theorem gives the following two equations 



and 

K-i 



[(1 - „) + (u - 0(1 - r-)]"" = E (^', 0 [(" - • (29) 

Substituting Equations 28 and 29 into Equation 27 gives 

• / (r*)'-+'*~^(l-r*)'+^-^tfr*. (30) 
Using the definition of the Beta /unction, the integral in Ecjuation 30 can be written as 



Consequently, Equation 30 can be written as 
rfY V V r(a + r(a -f r)r(/j + 

• (f ) Q (^'7 *) ^'"'(^ - «)'^~'"'(» - ly-^'' ■ (31) 

The gamma functions on the right side of Equation 31 can be written lis 

r{a)r{/?) r(« + /? + r + .^) (Q + ;^)r+/ 

where (x)j = + 1) . ■ . (a- + j - 1) and {x)q = 1. Substituting Equation 32 hi Equation 
31 gives 

p^(-v = .) = (f ) E E (:) ('\: ') (S^'-o - - • 



r=0 ji=0 

10 
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Let r' ~ i - r and s' ~ K ~ i — s. Changing the index of the first summation in 
Equation 33 from r to r' and the index of the jecond summation from a to s' produces 

i h'-i 



Px[X = i) = 



=<-.'EE(:o(;)f;') 

The terms in Equation 34 involving combinations can be written as 



T r r 



1 -u 
u-l 



(34) 



K\n\(K-i\ K\ t! (A' - i)\ 

i)\r')\ 5' ) i\{K-i)\ rn{i-r')\ sn{K -i-s')\ 




1 



r'!^'! {i'-r')\{K -i-s')\ 
[K ~ry.{K~s')\ 



r'\{K-T')\ s'\{K-s')\ {i-r')\{K -i-s')\ 
{K~r')\{K~sJ. (K\ (K 



(36) 



Substituting Equation 35 into 34 and rearranging the terms gives 

(cv)i-r' 



Pr(A' = e) = {u - If ^ 



{i - r')l 



r'=0 s'=0 

(K -r')l{K-s')l 



0 (r^)' 

CO (H) 



{K - i - s') 



(36) 



A'!(a + ^)K-r'_,' 

Lord (19C4) suggests using Equation 36 to calculate the observed score distribution when 
k = 0. Equation 36 can be calculated as the diagonal tenns of a matrix product of five 
(A' -}- 1) X (A' + 1) matrices (A/,, A/2, • . . , A/5) with the matrices based on the bracketed 
terms in Eqiiation 36. The lower triangular matrix Mi is given by 

0 0 ... 0 \ 

. 0 
. 0 



A/, = 



' 0? 
1! 



2! 



1 <'>>y 



(«)o 

0! 
(a)i 

1! 



0 
0 



0! 



i<*)K -I {oi)k-i 
(A'-J)! (A-2)? 



0! / 
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The diagonal matrix M2 is given by 



0 

A/2 = 0 
0 

The matrix M3 is given by 
/ 



0 

0 



M3 = 





Jv!(K-l)! 






(A'-l)!ff! 


(/v'-DU/f-l) 




• 




• 

2!(/C-l)! 


K!(a-»-/J)3 




I'K"! 










0 







0 

0 



A!(/t-2)! 



2!(A-2)! 
0 



The diagonal matrix M4 is given by 



M4 = 



0 
0 

0 



0 

0 

* 

0 



0 
0 



0 



The matrix A/5 is given by 



Ms- 



(A-l)! { A-2)! 



0 



i2k 

0! 

0 
0 



0 
0 
0 



\ 



A-?0! V 
0 



0 
0 

0 



/ 



\ 



0 
0 
0 



01 ' 
0 

0 
0 

0 / 
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The elements of each of the matrices Mi can be computed by first calculating an easily 
computed initial element and computing the o^her elements as a simple factor times an 
adjacent element. In addition, not all elements of the matrices need to be se^>arateh' 
calculatetl. For example, for matrix M\ only the first coluim needs to be computed 
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Four-parameter Beta Compound Binomial 

because the numbers in the other columns are a subset of the numbers in the first column. 
The fitted probabiUty for raw score i is the {i -f l)-th diagonal element of the matrix 
product M1M2M3M4M5 times {u — /)*'. 

Case 2: k > O.W hen k > 0, substituting the conditional error distribution given in Equa- 
tion 3 into Equation 1 gives the observed score distribution as 

Pr(.Y = i)= po{i) - J\t{1- t) [p{i j A' - 2, r) - 2p(i - 1 1 A - 2, r) 

+p{i - 2 I A - 2, r )] s(r I a, fi, /, u)dT , (37) 

where po(0 is the probability given in Equation 25 of raw score i when the conditional 
error distribution is binomial, p{i | n, r) is a binomial probabiUty as defined in Equation 
3, and g{T | a, j3, 1, u) is the four-parameter beta true score distribution given by Equation 
2. The terms r(l — r)p(j ) A' — 2, r) in Equation 37 can be written as 



^ 0 + l)(A-l-i) A! , 

A(A-l) (> + i)!(A-l-j)!'^ 

Substituting the result in Equation 38 into 37 produces 

Pr(.Y = i)= Poii) - j^^j^^^^ [{i + 1) (A - 1 - i)p{i + 1 I A, r) 

-2i{A - i)p{i I A, r)^{i-\){K + 1- i)p{i - 1 I A,r)] g{T | aj, l,u)dr. (39) 

Define p2(i) as 

pl{j)=j{K~j)poU)- (40) 
Using Equation 40, Ecjuation 39 can be written as 

Pr(A- = /) = p,{i) - j^-A—-[p'^{i - 1) - 2pUi) + pW + 1)] - (41) 

13 
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Four-parameter Beta Compoimd Binomial 

The procedure for calculating the observed score distribution for a value of ^ > 0, 
say ko, is to first calculate the observed score distribution assuming Ar = 0 using Equation 
36 (producing po{i),i — 0, . . . ^K) and then to use Equation 41 to calculate the observed 
score distribution for k = ko. 

When k > 0 some of the probabilities computed using Equation 41 may be negative. 
When negative values are computed they are usually very small in magnitude. 

Classification Consistency Indexes 

In this section the results presented previously will be used to obtain formulas for 
computing the classification consistency indexes described by Hanson and Brennan (1990) 
assuming the four-parameter beta compound binomial model holds for the test score in 
question. The two types of classification consistency indexes discussed by Hanson and 
Brennan (1990) are considered separately. Fir«t, the calculation of the probability of a 
consistent decision (coefficient p) and coefiicient k is discussed. Next, the calculation of 
the false negative and false positive error rates is discussed. 

Probability of a Consistent Decision and Coefficient Kappa 

The probability of a consistent decision (coefiicient p) is defined in terms of the bivari- 
ate distribution of scores on two independent administrations of a test. Assuming the test 
score random variables on two independent adiiunistrations of a test (denoted and X2 ) 
follow the four-parameter beta compound binomial model the probability of a randomly 
chosen examinee obtaining raw scores i and j on the two test administrations is (using the 
notation of Equation 1) 

Ft{Xi =i,X2=j) = Pr(X, =i\t, k) ?t{X2 = j | r, k) g{T \ «, /?, /, u) dr . (42) 
The classification consistency index p is c< ined as 

14 
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Four-parameter Beta Compound Binomial 

where Xo is the raw score cut point. 
Coefficient k is defined as 

^ _ P-Pc 

- = 13^^ (44) 
where is the probability of consistent classification by chance given by 

The bi variate distribution of scores on two independent administrations of a test given 
by Equation 42 can be computed based on the conditional error and true score distributions 
estimated from a single test administration. Thus, the first step in computing the bivariate 
distribution given in Equation 42 is to estimate the parameters of the four-parameter beta 
true score distribution using data from a single test administration as previously described. 
After the parameters of the true score distribution have been estimated, the bivariate 
distribution given in Equation 42 can be computed. 

Case 1: k^O, For computing the integral in Equation 42 the case is first considered in 
which the conditional error distributions [Fr(.Yj = i | r,Ar) and Pr(X2 = j | r,Jl-)] are 
binomial (A: = 0). In this case Equation 42 can be written as 

Pr( JCi = ^ X2 = i ) - ^ rr{X, ^z\r,k = 0) VriX^ = j | r, A^ = 0) gir j a, 0, /, «) dr 

[i^j)'- 'i^-''f''~'~'9ir\a,l3J,u)dT (45) 
The integral on the last line of Equation 45 is the probability of a test score of on a test 
with 2K items where the test score distribution is assumed to follow the four-parameter 
beta binomial model with the same true score distribution as the test .scores X, and X2. 
This integral can be calculated using Equation 36 (using the estimated parameters of the 
trur score distribution comj)utt'(l from a singh- test administration). 

2^ 
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Case 2: A; > 0, When Jt > 0, substituting the conditional error distribution given in Equa- 
tion 3 into Equation 42 gives the bivaxiate distribution of scores on two independent 
administrations of a test as 

Pr(-Yi =i,X2 = j) - J^'{p{i 1 K,T) 

~kT{l-T)[p{i\K-%T)-2p{i-l\K-2,T)+pii--2\K-2,r)]] 

• {pU I K,r) - fcr(l - T)[p{j \K-2,t)- 2p{j - 1 | A' -2,r) +p(j -2\K~2,t)]] 

-g{T\aJ,hu)dT, (46) 

/where p{i \ n, r) is a binomial probability as defined in Equation 3. Using the equality 

r(l - r)p(j I K - 2,r) = (l+^ffizii pO' + 1 1 K,r) 
from Equation 38 and expanding the terms in Equation 46 gives 

PviXi=i,X2 =j) - ^ I A',r)p(j I A,r) 

+ 1) [-(^" + 1) - i - + 1 1 a; r)p(j I A', r) 

-(i - 1)(A - i + - 1 1 A',r)p(j I K,t) 
+ 1)(A - j - l)p{i I A',r)p(j + 1 ! A',r) 
-{;■ - 1) (A - J + l)p{i I K%T)p{j - 1 I AT, r) 
+2 [(A - 0 / + (A - pii 1 A, r)p(; | A, r)] 

'^ A^(K-l)^ [(^ + 1)0 + 1) (A' - t - 1) (A - j - l)p(i + 1 I A, T)p{j + 1 I A, r) 
+(^ - 1)0' + 1)(A - t + 1){A -i - l)p{i ~ I I A,r)p(i + 1 I K,t) 
+(i + l)(j - 1)(A - i - 1)(A -i 4- l)p(i + 1 I A,r)p(i - 1 I A,r) 
- l)(i - 1)(A - i + 1)(A + l)p{i - 1 I K,T)p{j - 1 I A,r) 
-2 i (j + 1) (A - i) (A - j - l)p{i I 7v, r)p(i + 1 1 A, r) 
-2 f { j - 1) (A - 0 (A - j + l)p(e I K, T)p{j - 1 I A, r) 
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-2 (I + 1) j ( A' ~ i - 1) (A' - )p{i + 1 I A', r)p{j I A', r ) 
-2{i - l)j {K - i + l){K-j)p{i ~ 1 j A',r)p{i I A',r) 

+4 i j {K -i){K-j) p( I I K\ r ) p{ j J A', r )] | r 1 a, /?, /, tx) t/r , (47) 

Using Equation 45, the integral in Equation 47 can be evaluated to giw 

Pr(-Y, = f,.Y2=i)=po(f,;) 

'^ K{K-l) { r<{K-l^ ^' ~ - 1) (A' - I + 1) (A' - i + 1)] poii - XJ - 1) 



-f(i-l){A'-i + l) 
k 



-2k 



j(A'-j)-lJ Poii- I J) 



(A'-DA- 

[(i - 1)0' + 1)(A' - i + l){K-j - l)]po(i - l,i -f 1) 
-2k "J 



A-(A'-l) 

+0-i)(/i:-i + i) 



(A- 



+ {2 [(A- - ^)i^iK-J)J] ^ [4,-,-(A- - i)^K -i)]}po(r, j) 

( A^-l)A' '^^ ~'^"'^J + 



+(j + l)(A'-i-l) 
A- 



[(^" + - 1) (A' - / - 1)(A' + 1)] po(f + 1,; - 1) 
-2k 



I<{K~l) 

+ (2 + l)(A'~/-l) 



Poii + 



(A'-l)A- 

-^ A'(A--r) ^' ~^~l)iK-j- l)po{t + l,i + 1)1 , (48) 

where Poiij) is the probability of scores i and j on two independent administrations of 
the test when the conditional ernw distribution is binomial (given by Equation 45). 

When A' > 0 Equation 45 is used to first calculate the values of po(i, j), i, j = 0, . . . , A' 
which are used in Equation 48 to calculate tlie bivariate distribution of scores on two 
independent administrations of the test from which coefficient p and coefficient k can be 
computed. 
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Four-parameter Beta Compound Binomial 

False Positive and False Negative Error Rates 

The false positive error rate is given by 



V r?TiX = i\T,k)giT\a,^J,u)dT, 



m 



where Tq is the true score cutoff. The false negative error ratr is given by 

^ r Pt{X = i\ r,k)g{T\a,0J,u)dr , (50) 

isQ •''TO 

Case 1: k = 0. The case is first considered in which the conditional error distribution is 
binomial (it = 0). Based on a series of steps analogous to those which produced Equation 
30 from Equation 25 each term in the siun given in Equation 49 can be written as 

p Pr(X = i I T, Jfe = 0) ff(r j a, ^, /, u) dr 

=Jfe EE or: 

, / (^*)r+a-l(i_^*).+^~lj^*^ (51) 

where r* is defined as in Equation 22 and Tq ~ (tq - l)/[u ~ I). Using the definition of the 
incomplete beta function {/j.(a,6)] 



/.(«,ft)= r t"-' {I - tf-' dt , 

£f\a,b} Jo 

Equation 51 can be written as (analogous to Equation 31) 

Pr(X = ,\T,k = 0)s(r I a,/),/,«)dr = g g ^^^-^ _____ 

■I.;(a + r,/? + s) ) Q r ; ^) /-'{I - «)^— '(« - /)•^^ (52) 

In a manner analogous to that by which Equation 36 was obtaineci from Equation 31, 
Equation 52 can be rewritten as 

r Vr{X = i\r,k=0)giT\ a, /, u)dT = (« - /)^' L 
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Four- parameter Beta Compound Binomial 



A'!(a-f 



C') («-/) 



(53) 



Equation 53 can be used to calculate the terms in the sum of Equation 49 to give the false 
positive rate when k = 0. 

Case 2: k > 0. When A' > 0 an adjustment formula analogous to Equation 41 can be used. 
If Zi is the value of the integral in Equation 53 when k — 0 then the value of the integral 
when A- > 0 is 



k 



" /v(A'-l) 



(£- - 1) ( A' - i + ] ) r,-, - 2iiK - i) Zi + {i + 1) ( A' - t - 1) 



■ (54) 



The false negative error rate can be calculated using the values given by Equation 63 
(or when k > 0 Equation 54) since 

Ji J To 



SO that 



/ Pr(.Y ^ i I r,k)g{r \ oJi,lu)dT = 1 - / 'Pr(A' = i \ r,k)giT | a, pf,/, ti)t/r , 

J To J I 



(55) 



Values of the false positive and false negative error rates are computed based on the 
values from Eciuation 53 (or Equation 54 when k > 0) and 55, respectively. 

Computer Programs 

A set of functions written in the C language implementing the procedures presented 
in this paper for producing method of moments estimates of the true score distribution 
and observed score distribution under the f[>ur- parameter beta compound binomial model 
and computing classification consistency intlexes are available from the author. These 
functions can be compiled on any computer witli an ANSI C compiler. 
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Four-parameter Beta Compound Binomial 

Macintosh programs which use these routines are also av^lable from the author. One 
program calculates the method of moments Mtimates of the four-parmneter beta compoxmd 
binomial model, displays the results in a text window and produces graphic displays of the 
fitted distributions. This program also computes fitted raw score distributions based on a 
log-linear model (Holland and Thayer, 19S7). Another Macintosh program computes and 
displays the classification consistency indexes described for three models: the beta bino- 
mial, the four-parameter beta binomial, and the four-parameter beta compound binomial 
This program also displays the fitted models graphically. 
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